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ABSTRACT 

Motivated by the co mments of IGoodmaiil (j2009( ) on our paper concerning thermal tides 
(j Arras fc Socratejl2009af ). we have studied an idealized problem to under stand the global respons e 
of a completely fluid gas giant planet to thermal forcing at the surface (|Arras fc Socrates! l2009bh . 
Our findings disagree with the main claims in Goodman (2009). We find that significant quadrupole 
moments can indeed be induced as a result of thermal forcing. Furthermore, we find that it is possible 
for the orientation of the quadrupoles to be such that the planet is torqued away from synchronous 
rotatio n. Given these results, we believe our proposed thermal tide mechanism (lArras fc Socratesl 
l2009a| ) provides a viable scenario for generating steady-state asychronous rotation, inflated radii and 
possibly eccentric orbits of the hot Jupiters. 

Subject headings: planets 
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1. THE ISSUE 

Gold & Soter (1969; GS from here on) originally 
developed the idea of thermal tide torques to explain 
Venus' asynchronous spi n rate. Drawing on their work, 
lArras fc Socratei ()2009a[ ) assesed the importance of ther- 
mal tides for the hot Jupiters. Using the simple GS pre- 
scription for the quadrupole moment, they found that 
thermal tides could induce large asynchronous spin, and 
generate tidal heatin g rates more than sufficient to power 
the observed radii. iGoodmanI ()2009[ ) correctly pointed 
out that the GS ansatz does not faithfully represent the 
fluid motion induced by timc-dcpcndent heating in a 
completely fluid atmosphere. He argued that the induced 
quadrupole moment would be many orders of magnitude 
smaller than the GS value, and with an orientation that 
would act to synchronize the spin, opposite th e GS r esult. 

Motivated by the criticism of IGoodmanI (|2009| ). we 
attempted to carefully analyze a simplifled problem 
which captures the basic physics of thermal tide exci- 
tation in fluid planets. Our results are presented in 
(| Arras fc Socratej l2009bD . From here on, and unless 
stated otherwise, all references to equa t ions, fi gures and 
paper sections are to lArras fc Socratesl ()2009bf ) . 

In this note we compare our solutions to the fluid equa- 
tions (Ar ras fc Socrates 2 009b ) to the argumen ts pre- 
sented in .GoodmanI (|2009[ ). As iGoodmaiil (|2009D is un- 
published, we quote the text from his paper posted on th e 
Cornell University astro-ph archive (http://arxiv.org/). 
We then comment on the accuracy of the GS formula 
for the quadrupole moment. Contemp oraneous with 
the posting of Arras fc SP-Cjates (2009a), iGu fc Ogilvid 
(Ii009) published results on a related problem concerning 
thermal forcing of hot Jupiters. We briefly comment on 
the assumptions and results in these two papers. 

1.1. \Goodmari MM ) 
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The heart of Goodman's argument is contained in the 
fourth paragraph on his page 1: "A jovian planet, be- 
ing gaseous, lacks elastic strength. The excess column 
density of the colder parts of the atmosphere is counter- 
balanced to the degree that hydrostatic equilibrium holds 
by an indentation of the convective boundary and a redis- 
tribution of the cores mass toward the hotter longitudes. 
Insofar as the radial range over which mass redistribu- 
tion occurs is small compared to the planetary radius, 
the thermal tide therefore bears no net mass quadrupole. 
The torque on the atmosphere is opposed by a torque on 
the upper parts of the convection zone." 

Through these intuitive arguments, Goodman realized 
that the Gold and Soter approximation ignores the fol- 
lowing fact: though there is a flow from hot to cold at 
high altitudes, there is also a return flow at lower alti- 
tudes. In §5.1, we explicitly calculate the pattern of such 
a flow, by directly solving the fluid equations in the limit 
that inertia is ignored. The fact that this basic flow pat- 
tern is not contained in the derivation of the GS formula 
is a serious conceptual shortcoming. 

The analysis in §5.1 examines the limit of zero forcing 
frequency, and departs from Goodman's aforementioned 
paragraph in two respects. First, in the limit of zero 
forcing frequency, figure 3 shows that the return flow 
need not extend as deep as the convection zone. The 
bulk of the fluid motion is conflned in the radiative zone, 
near the photosphere of the starlight. Second, in §5.1 we 
do not find that density perturbations at high altitude 
are compensated by density perturbations of the opposite 
sign at lower altitude. Rather, we find that the density 
perturbation, and hence torque, is identically zero when 
fluid inertia is ignored. 

Given that density perturbations are zero in the limit 
of zero frequency, our next step was to derive flnite 
frequency corrections. Eq. 45 in §5.2 shows that fi- 
nite forcing frequency corrections give rise to a nonzero 
quadrupole moment. Have we violated hydrostatic bal- 
ance, assumed by Goodman, by including finite fre- 
quency? The answer rests on a technical detail, which 
may be important in future investigations. The equation 
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of hydrostatic balance dP/dz — —pg is obtained from the 
radial momentum equation by throwing away the inertial 
terms. Even if we were to throw away the inertial terms 
in the radial momentum equation, but kept them in the 
horizontal momentum equations, we would still find a 
nonzero quadrupole moment of the correct sign, although 
its magnitude would be slightly different (throwing away 
the —1 in the parenthesis in eq. 45 changes the prefactor 
from 4 to 5). 

By what factor should the GS quadrupole moment be 
reduced due to the "isostati c compensation" from the 
return flow? IGoodmaiil (|2009f ) argued above that the re- 
duction factor should be a power of H/R. By contrast, 
solution of the fluid equations in the limit of small forc- 
ing frequency, ignoring gravity waves, finds a frequency 
dependent reduction factor ~ 4(cr/A^)^ (see eq. 40). Al- 
lowing for gravity waves, the calculations in figures 4 and 
5 show that the response is larger than eq. 45 by 1-3 or- 
ders of magnitude in the relevant period range 1 day - 1 
month, due t o the excitation of gravity waves. 

iGoodmanI ()2009D clarifies the range of forcing period 
over which the quadrupole moment should be isostati- 
cally compensated on his page 2: "Whereas terrestrial 
isostasy operates on such long timescales that rock be- 
haves as fluid, the corresponding timescale for gaseous 
planets is dynamical, hence less than the tidal period." 

One of the key results of Arras & Socrates (2009b) is 
that low radial order gravity waves dominate the overlap 
with the thermal tide forcing. Hence the low frequency 
limit in eq. 45 docs not apply until forcing periods ~ 1 
month, comparable to or longer than the forcing peri- 
ods of interest (Arras & Socrates 2009a). Note that this 
surprising result is completely different from the case 
of gravitational forcing of incompressible fluid bodies, 
where the low frequency limit applies below the charac- 
teristic dynamical frequency {GMp/RpY^^ ^ (hours)"^ 
for a gas giant plane t. 

Lastly, IGoodmaiil (j2009l ) discusses the orientation of 
the induced quadrupole on his page 2: "Thus, the tidal 
torque claimed by Arras & Socrates (2009) vanishes to 
first order in the density variations of the thermal tide. 
To the next order, the quadrupole moment of the thermal 
tide aligns with the hottest and most distended parts 
of the atmosphere, because mass elements are weighted 
by the squares of their distances from the center. This 
will lead to a torque of the opposite sign to that of Ail, 
hence driving the planet toward synchronous rotation. 
Similarly, the phase lag of the thermal tide associated 
with an orbital eccentricity will affect the orbit only to 
second order, and will tend to circularize the orbit." 

The analytic solution to the fluid equations in the low 
frequency non-resonant limit (eq. 45) has the correct sign 
to drive asynchronous rotation, contrary to Goodman's 
claim. Including the effect of gravity waves, figures 4 
and 5 show that the sign of the quadrupole can alternate 
with forcing frequency. These sign changes are due to 
both the Lorentzian factors in eq. 60, as well as the 
signs of the quadrupole moments for individual modes 
(figure 6). Even in this more complicated case, frequency 
ranges still exist where the thermal and gravitational tide 
torques may oppose each other, leading to an equilibrium 
spin state. 

In summary, IGoodmanI ()2009[ ) correctly points out 



the deficiencies in the Gold and Soter approxima- 
tion employed by Arras & Socrates (2009a). How- 
ever, the solutions to th e fluid equations presented in 
lArras fc Socratei (|2009b[ ) differ both qualitatively and 
quantitatively from the basic picture outlined in his work. 
Consequently, we disagree with Goodman's criticism of 
Arras & Socrates (2009a) i.e., that thermal tides cannot 
lead to asynchronous spin and eccentric orbits. 

1.2. the Gold & Soter approximation, and the 
calculations of Arras & Socrates (2009a) 

Gold and Soter's ansatz involves a major assumption: 
that the fluid elements remain at roughly constant pres- 
sure, so that density perturbations are related to temper- 
ature perturbations by 5p/ p = —5T/T = — As/cp. From 
eq. 29, we see this is indeed true if the 5p and terms 
can be ignored. For low frequency forcing, ignoring the 
5p term may be a good approximation, but for fluid at- 
mospheres we have seen it is not a good approximation 
to ignore the term. If there is a solid surface, and the 
boundary condition at this surface is = 0, then the 
Gold and Soter ansatz may, in fact, hold. This might be 
realized if the heat was all deposited at or near the solid 
surface, rather than well above the surface. 

On a more practical level, for fluid planets with surface 
radiative layer, we have found the Gold and Soter ap- 
proximation overestimates the quadrupole moment and 
torque by more than an order of magnitude for the cal- 
culations in this paper. Since the steady state planetary 
radii powered by tidal heating found by Arras & Socrates 
(2009a) were rather large compared to observed planets, 
the reduction in torque found in this paper may bring 
their theory into better agreement with observations. 

To summarize, we have found that the Gold and Soter 
ansatz for the quadrupole moment is qualitatively, but 
not quantitatively correct. It may be viewed as a conve- 
nient order of magnitude estimate. 

1.3. \Gu & Oailvi2 MM ) 

iGu k, Qgilvid (|2009!) studied perturbations to a hot 
Jupiter atmosphere induced by time-dependent radiative 
heating due to asynchronous rotation. Their primary 
result is that waves excited by the thermal forcing can 
radiatively damp and transfer angular momentum verti- 
c ally in the atmospher e . 

iGu fc Qgilvid ()2009D did not study if net quadrupoles 
could be induced in the atmosphere. By working in 
plane parallel geometry and requiring the pressure per- 
turbation to vanish at the base of the grid, they set the 
quadrupole to zero by hand. As the goal of their paper 
was to study differential rotation induced by the damp- 
ing of downward propagating gravity waves, it is likely a 
good approximation to work in plane parallel geometry, 
ignoring the possible existence of net quadrupoles. Note 
that Arras & Socrates (2009b) focus on the complemen- 
tary issue of net quadrupoles, ignoring possible driving 
of differential rotation. 

2. SUMMARY 

In summary, the results in this paper confirm that 
quadrupole moments of the correct sign and approx- 
imately correct magnitude may be induced by time- 
dependent insolation, confirming the basic assumptions 
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of Arras & Socrates (2009a). A future study will use the 
results from this paper in concert with a thermal evolu- 
tion code for the hot Jupiters (Arras & Bildsten 2006; Ar- 
ras & Socrates 2009a) to construct more detailed steady 
state solutions for the planetary rotation and radius, and 
the orbital eccentricity. 
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